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Preface 

This  thes-.b  represents  my  attempts  to  determine  the  differ¬ 
ential  game  barrier  for  two  spacecraft  in  near-earth  orbit.  The 
vehicle  parameters  are  in  general  viewed  in  two  dimensions  while  the 
barrier,  due  to  the  dimensionality  of  the  linear  state  equations,  is 
constructed  in  four  dimensions.  Once  determined,  the  barrier  is 
analyzed  to  determine  whether  it  divides  the  playing  space  into  two 
separate  zones:  an  escape  zone  and  a  capture  zone. 

I  would  like  to  express  my  sincere  appreciation  to  the  faculty 
of  the  Air  Force  Institute  of  Technology  for  their  assistance  in  the 
research  and  preparation  of  this  thesis.  I  would  especially  like  to 
thank  my  advisor,  Major  Gerald  M.  Anderson,  for  his  guidance  and 
suggestions  toward  the  completion  of  this  study  and  also  Professor 
Edstrom  of  the  Mathematics  Department  for  devoting  his  time  to  help¬ 
ing  solve  many  differential  equation  problems.  I  also  want  to  express 
my  heartfelt  gratitude  to  my  -wife,  Marsha,  for  her  undying  devotion, 
unfailing  enthusiasm,  and  unceasing  encouragement.  Also,  I  want  to 
thank  the  Lord  for  this  opportunity  He  has  provided  and  for  His  inval¬ 
uable  guidance,  insight,  and  motivation  without  which  this  study  would 
not  have  been  completed. 


Victor  W.  Grazier 
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I*  Introduction 

Background 

The  advent  of  the  space  age  has  thrust  man  into  orbital  space- 
craft  and  outward  to  the  moon.  This  intense  activity  in  space  has 
made  it  seem  inevitable  that  man  would  arm  his  space  vehicles  in  an 
effort  to  gain  power  or  influence  on  earth.  In  an  effort  to  preclude 
such  actions,  many  nations  joined  together  in  signing  the  Treaty  of 
1967,  an  agreement  to  keep  space  neutral.  Although  man's  intentions 
were  commendable,  the  past  has  shown  dramatically  that  dreams  of 
power  can  overcome  any  barriers,  and  unfortunately,  the  barrier  to 
arms  in  space  must  be  viewed  in  this  light.  While  vigorously  support¬ 
ing  the  efforts  to  keep  space  neutral,  the  United  States  must  continue 
its  research  in  the  area  of  orbital  defenses  in  the  eventuality  that  some 
aggressive  government  attempts  to  gain  control  of  the  space  around  the 
earth. 

Problem 

This  study  investigates  the  region  surrounding  two  orbital 
spacecraft  engaged  in  a  pur  suit -e%rasion  game,  in  this  game,  the 
pursuer  desires  capture  while  the  evader  desires  to  escape.  Termi¬ 
nation,  then,  if  it  is  possible,  results  when  the  pursuer  forces  the 
evader  to  come  within  a  specified  distance  of  him.  This  distance  is 
the  radius  of  the  terminal  surface  which  is  centered  on  the  pursuer. 
From  this  surface,  trajectories  are  generated  backward  in  time  in 
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order  to  provide  information  concerning  another  surface  known  as 
the  "barrier.  " 

The  barrier  is  a  surface  which  results  from  the  application  of 
differential  game  theory.  This  surface  generally  divides  the  playing 
space  of  interest  into  two  zones:  a  capture  zone  and  an  escape  zone. 
Starting  points  from  which  capture  'an  not  occur  lie  in  the  escape 
zcne.  Those  from  which  escape  is  not  possible  lie  in  the  capture 
zone.  Starting  points  on  the  barrier  itself  will  remain  on  the  barrier 
for  a  neutral  outcome  as  long  as  both  players  use  their  optimal  strat¬ 
egies. 

Current  Knowledge 

The  use  of  differential  game  theory  in  numerous  fields  of 
endeavor  began  in  earnest  in  the  late  1950's.  A  major  contributor  to 
this  theory  was  Rufus  Isaacs  (Ref  1).  Although  differential  game 
theory  has  been  investigated  extensively  over  the  past  decade  for 
applications  to  air-to-air  engagements,  relatively  little  work  has  been 
done  on  the  use  of  this  theory  for  orbiting  spacecraft,  R.  E.  Wcng 
(R.ef  Z )  used  differential  game  techniques  to  investigate  fixed  time, 
minimum  final  separation;  assuming  that  the  encounter  took  place  in  a 
constant  gravity  field.  In  his  master's  thesis,  R.  H.  Woodward  (Ref  3) 
investigated  the  free  time,  minimum  final  distance  problem  in  an 
inverse  square  gravity  field.  To  date,  no  one  has  investigated  differ¬ 
ential  game  barriers  for  two  spacecraft  in  orbit. 
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Scope 

The  construction  of  the  barrier  requires  numerous  assumptions 
(some  of  which  are  listed  in  the  next  section)  along  with  many  approxi¬ 
mations.  These  are  necessary  to  reduce  the  dimensionality  of  the 
problem  to  something  that  can  be  handled,  while  also  removing  com¬ 
plexities  that  do  not  greatly  affect  the  outcome  of  the  problem.  In  this 
manner,  the  closed-form  solutions  may  be  found  which  ultimately 
will  allow  analysis  of  the  barrier  itself.  Once  the  barrier  has  been 
constructed,  the  objective  is  to  analyze  this  barrier  to  define  the 
conditions  necessary  for  establishing  whexu  ea.ch  starting  point  is 
located.  In  other  words,  by  applying  these  conditions,  it  should  be 
possible  to  determine  whether  given  starting  points  are  in  the  escape 
or  the  capture  zone.  Surely  other  objectives  could  be  sought;  however, 
due  to  time  iimitatio  is,  this  study  will  concentrate  only  on  the  above 
search. 

Assumptions 

1.  Both  vehicles  thrust  continuously..  If  the  evader  stopped 
thrusting,  capture  couru  readily  occur.  If  ihe  pursuer  quit  thrusting, 
the  evader  could  rapidly  escape. 

2.  This  pursuit-evasion  game  takes  place  in  an  inverse- 
square  gravity  field, 

3.  Both  vehicles  have  constant  mass.  This  entails  the  further 
assumptions  that  the  time  period  for  the  pursuit-evasion  game  is 
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sufficiently  short  and  that  the  thrust  is  very  low.  Variable  mass, 
however,  would  be  an  area  for  further  investigation  in  future  studies. 

4.  Control  is  provided  by  varying  the  thrust  vector  direction. 
The  control  angle  itself  is  measured  from  each  vehicle's  local  hori¬ 
zontal. 

5.  Each  player  has  perfect  knowledge  of  the  position  and  capa¬ 
bilities  of  his  opponent.  Otherwise,  the  game  would  degenerate  with 
one  player  holding  an  obvious  advantage. 

Approach 

1.  Linearize  the  equations  of  motion  about  a  circular  reference 

orbit. 

2.  Normalize  both  the  linear  and  nor.-linear  equations  of  motion 
about  the  same  circular  reference  orbit.  Use  normalized  time  to 
develop  the  linear  and  non-linear  state  equations. 

3.  Apply  necessary  conditions  to  establish  boundary  conditions 
for  both  the  linear  and  non-linear  equations. 

4.  Compare  linear  trajectories  with  corresponding  non-linear 
trajectories  in  order  to  determine  the  validity  of  the  linear  approxi¬ 
mation. 

5.  Vary  parameters  around  the  terminal  surface  and  out  the 
trajectories  in  order  to  gain  information  on  the  barrier. 

6.  Develop  a  method  to  determine  whether  each  given  starting 
point  lies  in  the  escape  zone,  in  the  capture  zone,  or  on  the  bariier 
itself. 
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II.  General  Equation  Development 

The  objectives  of  this  chapter  are  first  to  develop  the  general 
equations  of  motion  for  a  spacecraft  in  orbit,  then  to  linearize  these 
equations  about  a  circular  reierence  orbit,  and  finally  to  normalize 
them  about  this  same  reference  orbit. 


Nonlinear  Equations  of  Motion 

Assuming  that  the  earth  is  spherical  and  then  using  a  polar 
coordinate  system,  the  equations  of  motion  for  a  spacecraft  in  orbit, 
as  foimd  in  any  basic  dynamics  textbook,  are  as  follows: 


JL  4  T  sin  Cf 
TZ  m 


(1) 


t$  +  Zrd  =  H-co.s  .  (2) 

m 

where  Ol  is  the  angle  between  the  thrust  vector  and  the  local  .horizontal 
as  in  Fig.  1. 

In  order  to  work  with  these  equations,  it  is  desirable  to  trans¬ 
form  them  into  expressions  concerning  the  tangential  and  radial 
velocities,  Xq  and  respectively.  This  can  be  accomplished  by  using 
the  following  relationships. 


(3a) 

(3b) 
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Fig.  1.  Orbital  Geometry 
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Through  differentiation,  expressions  for  r,  0,  and  may  be 
obtained  which,  when  substituted  into  Eqs.  (1)  and  (2),  yield  these 
equations  of  motion. 

V  2 

V  -  +  T  sin  ot 

r  r  rZ  m 


V 


V  v 

rtf 


T  cos  Qt 
m 


Linearized  Equations  of  Motion 

Linearizing  Eqs.  (3)  and  (4)  about  a  circular  reference  or'ci' 

yields 


hi  =  Av 
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A0  - 


V^Ar 

- T~ 

r 


as  well  as 


AV 


„.2VSAV0  ,VAr 

r  2 


2(1  Ar 


_  T  sin  Of 


m 


(5b) 


(6) 


«,V  .  4VrV0  .  Vr^ve  VrV9Ar  T  cos  a 

AV0  + - j - + - - - j - * - 

r 


(7) 


In  these  equations  the  terms  V  ,  V/»,  and  r  now  refer  to  the 

r  u 

reference  orbit.  The  delta  terms  refer  to  deviations  from  the  refer¬ 
ence  orbit  for  these  given  parameters. 

Assuming  that  the  reference  orbit  for  linearization  is  circular. 


V 

r 


=  0 


and 


(8a) 


or 

v*  •£ 


(8b) 


Substituting  these  expressions  into  Eqs.  (6)  and  (7)  gives  the  following 
linearized  equations  of  motion. 

Av  .*VY»  .V*  .TsinO.  ,, 

r  r  r2  m 


AVfl  =  -  AVrV0  +  T  cos  a 
°  r  m 


(10) 
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Normalizing  the  Variables 

At  this  point  it  is  desirable  to  normalize  the  parameters  about 
a  circular  orbit,  with  radius  tq  and  velocity  Vq,  so  that  significant 
figures  are  not  lost  due  to  large  differences  in  magnitude  between 
these  parameters.  The  velocity  magnitude  may  be  several  miles  per 
second  while  the  radii  may  be  thousands  of  miles. 

In  order  to  work  with  the  non-dimensional  parameters,  espec¬ 
ially  regarding  the  state  equation  development,  it  is  necessary  to 
introduce  a  non-dimensional  time  parameter  T  defined  by 


T  = 


r 


0 


t 


(ID 


By  means  of  the  chain  rule  for  differentiation. 


ill  =  -.h  !U 

d  T  dt  dr  V  dt 


(12) 


In  this  manner,  —  ^  -  V^.,  where  the  prime  represents  differentiation 
with  respect  to  T. 


Non-linear  Normalization 


Define  the  following  normalized  state  variables. 


X, 


=  JL 

(13a) 

r0 

V, 

-  — — 

(13  b) 

vo 

1 

Vn 

(13c) 
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X4  =  6  (13d) 

By  use  of  Eq.  (12),  the  following  normalized,  non-linear  equa¬ 
tions  of  motion  are  determined. 

Xj  =  X2  (14a) 


d(X2) _  r0 

[Vrl 

=  r° 

rv2 

9  _  P  +  T  Sin  Ci 

dT  V0 

v0j 

V  ^ 

0 

t  2  m 

.  r  J 

For  the  normalizing  orbit,  fi  =  V  ^  Tq,  so  that 


X_ 

Xl  =  — 
2  X1 


+  U  sin  fit 


(14b) 


X, 


p  ip 

where  U  =  —2 _ ,  a  constant  with  thrust  T  and  mass  m.  Differenti- 

Vm 

«.lion  cf  Eqs.  (13c)  and  (i3d)  with  respect  to  T  yield  the  remaining 
non-linear  state  equations. 


X->X, 

x'  =  -  -  ...  3  +  u  cos  a 
3  X, 


(14c) 


„_X3 


(14d) 


To  this  point  the  general  equation  development,  has  followed 
closely  the  work  of  R.  H.  Woodward  (Ref  3).  In  order  to  avoid 
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repetition  and  since  the  normalized,  linear  equations  involve  terms 
for  both  the  pursuer  and  the  evader,  this  development  will  be  included 
in  the  next  chapter  along  with  the  differential  game  barrier  theory. 
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TTT.  Differential  Game  Barrier  Formulation 

This  chapter  presents  the  theory  and  equations  of  the  differ¬ 
ential  game  barrier  formulation  for  a  pursuer  and  evader  engaged  in 
a  pursuit-evasion  game  in  orbit. 

Theory 

This  section  must,  of  necessity,  follow  closely  the  work  of 
others.  The  nomenclature  and  theory  of  both  Isaacs  (Ref  1)  and 
Bryson  and  Ho  (Ref  4)  are  combined,  with  the  resulting  development 
being  very  similar  to  that  used  by  P.  H.  Cawdery  (Ref  5), 

The  objective  of  this  development  is  to  find  the  minimax  solu¬ 
tion,  if  one  exists,  to  the  expression 

J  =  0[X(tf)]  (15) 

with  n-dimensional  state  vector  of  the  type 

X  =  f_[X,  u,  v,  t];X  (t0)  =  3So  (l6) 

and  terminal  constraints 

^[X(tf)]  =  0  (17) 

In  Eq.  (16)  u  and  v  represent  the  controls  associated  with  the 
pursuer,  P,  and  the  evader,  E,  respectively.  The  aim  is  to  determine 

zk  -i- 

the  pair,  u  and  v'1',  such  that 
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IVIW  \il\f 


J  (u’1',  v)  <  J  (u*,VV)£  J  (u,  V*) 


According  to  Isaacs,  the  saddle  point  solution,  J  (u*,  v"'),  is  the  Value 


of  the  game.  The  necessary  conditions  for  this  solution,  called  the 


solution  in  the  small,  are 


H(X,  X,  u,v,t)  =  \lf_ 


X  =  -  Hx  (X,  ^  u,  v,  t) 


where  H  is  defined  as  the  Hamiltonian  and  \  is  the  costate  vector. 

The  n- dimensional  costate  vector  X  is  subject  to  the  following 
transversality  conditions. 


X  (tf)  =  (tf )  +  v  &  (tf) 


H  (tf )  =  <f>t  (tf)  +  (t£) 


where  v  is  a  Lagrange  multiplier. 

In  order  to  determine  the  optimal  controls,  H,  of  Eq.  (20), 
must  be  minimized  with  respect  to  control  u  and  maximized  with 
respect  to  v  so  that 


H*  =  max  min  H  (X.  X.  u,  v,  t) 


Tn  order  for  the  solution  to  exist,  H  must  be  assumed  to  be 


separable.  In  this  case 
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max  min  IT  -  min  max  H  (24) 

V  U  U  V 

For  this  problem  H  is  assumed  to  be  separable.  This  means  that  H 
does  not  have  terms  containing  both  u  and  v  together.  This  insures 

j, 

a  saddle  point  solution  as  stated  in  Eq.  (18)  with  u  minimizing  and  v 
maximizing.  In  addition,  if  t  does  not  explicitly  appear  in  Eq.  (19)* 
then  H,  4>,  and  are  not  dependent  on  time  thereby  yielding  the  follow¬ 
ing  important  result  from  Eq.  (22). 

H(t)  =  H(tf)  =  0  (25) 

Games  of  Kind 

Isaacs  (Ref  1)  introduces  two  types  of  differential  games: 
games  of  degree  and  games  of  kind.  In  the  former,  which  may  be  con¬ 
tained  within  a  game  of  kind,  the  players  attempt  to  maximize  or  mini¬ 
mize  the  final  payoff  as  a  function  of  the  variables  of  the  game.  In  this 
type  of  game  the  payoff  is  not  determined  until  termination,  thus  pre¬ 
suming  that  the  game  does  actually  end.  In  games  of  kind  the  essential 
problem  is  the  achievement  of  termination  where  one  player 

desires  termination  while  the  other  does  uot. 

Termination  in  Isaacs'  game  of  kind  is  taken  to  mean  either 
that  the  evader  is  forced  to  come  within  a  certain  proximity  of  the  pur¬ 
suer  or  a  long  time  period  elapses  without  capture  occurring.  The 
proximity  has  been  formalized  in  the  playing  space  into  a  surface  C 
which  is  defined  as  the  terminal  surface.  For  capture  to  occur,  then, 
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the  evader  must  be  forced  not  only  to  this  surface  but  also  to  penetrate 
it.  Capture,  escape,  and  a  third  result,  that  of  a  neutral  outcome,  are 
explained  more  fully  below. 

In  this  thesis  the  pursuer  P  desires  termination  or  capture 
whereas  the  evader  E  desires  to  avoid  termination,  thus  to  escape,  or 
at  least  to  delay  capture  as  long  as  possible.  The  outcomes  of  the 
game  of  kind  then  fall  into  three  general  categories: 

a.  P  can  capture  for  all  starting  positions  in  the  playing  space. 

b.  E  can  escape  for  all  starting  positions. 

c.  Neither  P  can  capture  nor  E  can  escape  for  all  starting 
positions:  a  neutral  outcome. 

Category  c  is  the  one  of  interest  in  this  thesis  since  the  result 
demands  that  both  P  and  E  play  optimally.  If  P  were  to  momentarily 
play  nonoptimally,  F.  could  force  P  out  of  a  neutral  outcome  and  escape, 
likewise,  if  E  were  to  play  nonoptimally,  capture  would  result.  There¬ 
fore,  the  playing  space  is  divided  into  two  zones,  capture  and  escape, 
by  a  surface  which  contains  all  starting  points  for  which  the  outcome 
is  neutral.  This  surface,  then,  is  defined  as  the  barrier. 


The  Barrier  Development 

The  condition  that  must  be  met  for  a  barrier  to  exist  is  as 
follows : 


MC/GA/73A-3 


where  v  represents  a  vector  normal  to  the  terminal  surface.  This 
relationship  will  be  called  the  Barrier  Hamiltonian,  a  term  borrowed 
from  P.  H.  Cawdery  (Ref  5),  and  is  represented  by  H^.  Since  the  u 
vector  dotted  into  {_  equals  zero,  where  l_  equals  the  velocity  X,  the  two 
vectors  must  be  orthogonal.  Therefore,  at  that  point  on  the  terminal 
surface,  the  component  of  the  velocity  in  the  direction  of  is  zero, 
thus  penetration  can  not  occur.  Points  which  meet  this  criteria  have 
been  classified  by  Isaacs  as  the  boundary  of  the  usable  part  (BUP).  The 
usable  part,  where  Eq.  (26)  is  less  than  zero,  defines  the  portion  of  the 
terminal  surface  where  the  pursuer  P  can  force  termination.  Obvi¬ 
ously,  then,  the  nonusable  part  is  that  portion  of  the  terminal  surface 
where  Eq.  (26)  is  greater  than  zero.  At  these  points,  the  evader  E  can 
frustrate  termination  and  thus  escape;  therefore,  this  portion  is  non¬ 
usable. 

The  construction  of  the  barrier  involves  starting  with  the  initial 
conditions  at  the  terminal  surface  and  integrating  the  state  and  costate 
equations  backwards  in  time.  The  resulting  trajectories  provide  points 
required  for  constructing  the  barrier. 

Non-linear  Equation  Development 

The  non-linear  equations,  presented  earlier  as  Eqs.  (14),  are 
repeated  below  in  the  differential  game  format  for  the  evader.  The 
equations  are  identical  for  both  players  with  only  the  (E)  or  (P)  sub¬ 
scripts  to  duierentiate  between  the  two. 
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v 


X1E  =  X2E 


(27a) 


X 


k3E  1 


2E 


1  U  v  sin  d 


X 


IE  X1£‘ 


E 


(27b) 


X3E  =  “  UE  cos  Q'E 


(27c) 


IE 


X 


'  X3E 


4E 


(27d) 


‘IE 


A  like  set  of  equations  exists  for  the  pursuer  with  the  subscript  (E) 
replaced  by  (P). 

The  Hamiltonian,  then,  is  formed  by  adjoining  the  costates  tc 
the  state  equations. 


H  "  X1EX2E  +  A2E 


x3e‘ 


IE  X 


J _  +  UE  sin 

IE 


+  \ 


3E 


v  v 

"‘2EX3E 


X 


IE 


+  UE  cos  atE 


Xi  p 

+  A4E 

X1E 


'■  ^ipX'zp  +  A2P 


r  ? 

X3P 

~ -  ’ - W  +  Up  sin  «p 

AiP  XlpZ 


+  0 


X2PX3P  1  ,  X3p 

+  Up  cos  op|  4  X4P 


X 


IP 


‘IP 


(28) 


From  this  equation  the  barrier  costate  equations  may  be  determined  by 
applying  the  necessary  conditions  from  Eqs.  (20). 
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d  H 


dX 


IE 


X2EX3E2  .  2X2E  .  ^3EX2EX3E  +  >-4EX3E 

y  2  y  2  y.2  y  2 

X1E  X1E  X1E  X1E 


(29a) 


X1E  + 


A3EX3E 


X 


IE 


(29b) 


2A2EX3E  +  a3EX2E 
X1E  X1E 


X4E 

X1E 


(29c) 


X4E  ~  0  yielding  =  a  constant  (29d) 

Again,  a  like  set  of  equations  exists  for  the  pursuer  with  the  subscript 
(Pd)  replaced  by  (P). 

The  optimum  barrier  control  angles  for  the  pursuer's  and  the 
evader's  thrust  can  be  determ  :ed,  as  mentioned  previously,  by  mini¬ 
mizing  H  with  respect  to  u  and  maximizing  it  with  respect  to  v,  v/here 
Otp  represents  u  and  Ct  g  represents  v.  After  minimaxing  H,  the 
proper  signs  are  determined  by  applying  the  following  second  partial 
derivative  necessary  conditions. 

>  o  <  o  (30) 

dir  dv2 


The  maximizing  and  minimizing  of  H,  then,  results  in  the  following 


tan  a  £ 


=  X2E 
*3E 


(31) 
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and 


tan  Cl 


(32) 


By  applying  the  conditions  in  Eqs.  (30),  the  optimum  barrier  control 
angles  are  determined  to  be 


*3E 


cos  dr  = 


X2 

+  A3E 


and 


-  X 


2  P 


liiP- 


sin  dp  =  x2 

^A2P  A3P 


COS  01-1 


V 


'X2P  +  X3P 


(33) 


(34) 


These  conditions  may  now  be  substituted  into  the  state  and  costate 
equations  to  yield  the  TPBVP  as  functions  of  only  the  states  and  the 
costates  themselves. 


Linear  Equation  Development 

As  mentioned  in  Chapter  II,  the  linear  equation  development 
involves  both  the  pursuer  and  the  evader  almost  from  the  beginning  so 
the  entire  development  was  delayed  until  this  time.  Again  there  is  the 
requirement  to  define  state  variables,  re  being  the  linearized,  nor¬ 
malized  state  variables  listed  below. 


Ax 


IE 


AX 


2E 


Ar 


E 


AV 


rE 


(35a) 

(35b) 
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AX 


3E  7 


(35c) 


AX4E  =  A0 


(35d) 


A  like  set  of  terms  exists  for  the  pursuer  with  the  subscript  (E) 
replaced  by  (P). 

In  order  to  reduce  the  number  of  dimensions  in  this  problem, 
it  was  necessary  to  combine  the  pursuer's  and  the  evader's  state  vari¬ 
ables  into  one  difference  relationship,  now  represented  by  the  Y 
vector.  Fig.  2  depicts  the  linearized  parameters  which  have  been 
combined  into  Yj  and  Y4. 

The  linearized,  normalized  state  variables  then  become 


Yi“ 


^X1E  ~  ^X1P  =  ~  ^rp 


(36a) 


Y3  ^ 

Y4  = 


rE  -  AvrP 


-AV0P 


A0e  -  A0p 


(36b) 


(36c) 


(36d) 


Differentiating  these  pmiafions  with  respect  to  T  bv  use  of 
Eq.  (12)  and  replacing  the  dotted  terms  by  use  of  Eqs.  (3)  and  (4)  result 
in  the  following  linearized  state  equations. 

IS) 


Earth's  Center 
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rP 


AV, 


6E 


IV 

'0 


+ 


+  ro 

TEsin  aE 

Tpsin  op 

r  2  r  2 

l  E  P 

v0 

mE 

i m 

mp 

VAVrP 


Tecos  oe 


mE 


* 

Tpcos  ap 
mP 


ro 

r^E 

iv«pl 

r0 

V0E  ArE 

vo 

rE 

tm 

"P  _ 

'vo 

2 

L  rE 

rp2  J 

With  the  reference  orbit  identical  to  the  normalizing  orbit, 


r  = 


r 


0 


and 


'  n 
V 


=  V 


Using  these  expressions  and  by  once  again  letting  U  = 


the 


(37a) 


(37b) 


(37c) 


(37d) 


above 
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equations  reduce  to  the  following  state  equations. 

Y[  =  Y2  (38a) 

Y2  =  2Y3  +  Yl  +  UE  sin  -  Up  sin  0tp  (38b) 

Y3  =  -  Y2  +  Up.  cos  Clp  -  Upcos  Op  (38c) 

Yi  =  Y3-  Yj  (33d) 


These  equations  are  new  linear  in  state  only,  still  being  non-linear  in 
the  control  angle  O.  They  will,  however,  continue  to  be  referred  to  as 
linear  equations. 

This  development  for  the  linear  state  equations  again  is  similar 
to  the  work  of  R.  H.  Woodward  (Ref.  3). 

The  Hamiltonian  for  the  above  equations  is  again  formed  by 
adjoining  the  costates  to  these  state  equations. 

H  =  A  Aj  Y?  +  A  A2  £2 Y-  +  Y p  +  Up  sin  a,-.  -  UpsinOtpj 

+  A  A3  f-  Y3  +  Up  cos  Ctp  -  Upcos  otpl  +  AA4  [y3  -  Y jj  (39) 


By  once  again  applying  the  necessary  conditions  of  Eq.  (20),  the  follow¬ 
ing  costate  equations  are  formed. 

AA'  =  -  -2H-  =  -  A  A  +  A  A  (40a) 

1  d Yj  “  4 

A A^  =  -  AA1  +  A A3  (40b) 


22 


MC/GA/T3A-3 


A\'3  =  -  2A\2  -AX4  (40c) 

AX4  =  0  yielding  AX^  =  a  constant  (40d) 

The  optimum  thrust  angles  are  found  in  the  same  manner  a3  in 
the  non-linear  development  yielding 

AX, 

tan  a  =  —L  (41) 

E  Ax3 

and 

tan  Ct0  =  _ (42) 

P  AAj 

By  applying  the  second  partial  derivative  relationships  of  Eqs.  (30),  the 
following  optimum  barrier  control  angles  are  determined. 


sin  ar  =  +  A  X2 


V 


rAX?  +  AX3 


sin  a. 


(43) 


cos  a  -  +  ^  ^3 

ili 


y^T 


AX? 


=  cos  Oli 


(44) 


These  equations  demonstrate  that  when  both  the  pursuer  and 
the  evader  play  their  optimal  strategies,  their  thrust  angles  will  be 
the  same.  On  substituting  these  expressions  into  the  state  and  costate 
equations,  these  equations  reduce  to  a  TPBVP  again  as  functions  of 
only  the  states  and  costates  themselves. 
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Terminal  Conditions 

Non-linear  Development.  A  circular  terminal  surface  is 
assumed,  lying  in  the  plane  containing  the  pursuer,  the  evader,  and 
the  center  of  the  earth.  This  surface,  then,  provides  the  end  condi¬ 
tions  for  certain  state  and  costate  equations.  These  conditions  are 
determined  from  the  following  relationship  representing  the  assumed 
terminal  surface. 


where 

-  Xfp  =  R  sin  <f>  (46a) 

and 

X4E  -  X4p  =  R  cos  <f>  (46b) 

In  this  expression  is  defined  as  the  angle  around  the  terminal  sur¬ 
face  measured  counterclockwise  from  the  pursuer's  local  horizontal. 
The  terminal  surface  and  parameters  are  depicted  in  Fig.  3. 

The  vector  v_  mentioned  previously  as  normal  to  the  terminal 
surface  is  found  by  taking  the  gradient  of  Eq.  (45).  This  operation 
yields  the  following  terminal  state  and  costate  relationships. 
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V  [Eq.  (45)]  = 

rx1E 

“  X1P 

x 

sin  ip 

R 

aie 

X4E 

“  X4P 

*4E 

cos  <p 

R 

X1P 

-  X1E 

X 

-  sin  <p 

R 

A1P 

x4P 

"  X4E 

Kr, 

-  cos  <p 

- 

R  _ 

4P 

».  1 

as  well  as  A2E,  *3E>  ^2P»  ^3P  =  0  (47) 


Here,  vector  v  is  represented  by  Ajp.,  • .  . ,  A4E  A  jp>  •  •  •  *  A4p- 
The  remainder  of  the  terminal  conditions  may  be  found  by  using  the 
condition  tor  the  boundary  of  the  usable  part,  Eq.  (26), 


HB^*x  =  x1Ex2E  +  A2E 


rx  2 

A3E  i 


X1E  X1E 


4~  sin  Qp-> 


+  A 


3E 


X2£X3E  "|  A4EX3E 

+  UE  cos  dEi  +  +  ^\pX2p 


X1E 


IE 


+  A 


2P 


rx  2 

A3P 

i 

+  Up  sin  flip 

+  A 

X2pX3P 

-  - .  4-  TT  rrtS  a 

V 

v2 

3  P 

t  Up  Up 

L^ip 

AlP 

'J 

1 

■"  1  o 

L  -1 

+  ^4PX3p 

x1P 


(48) 
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This  equation  is  identical  to  the  Hamiltonian  of  Eq.  (28),  Therefore, 
HgCtf)  =  0  as  in  Eq.  (25),  and  using  the  terminal  conditions  in  Eq.  (47), 
this  reduces  to 


! 


X2£  "  X2P! 


sin  0  + 


X3E  X3P 


x1E  X1P 


cos  <j>  -  0 


(49) 


To  determine  the  thrust  angle  values  at  the  terminal  surface, 
tan  a  is  evaluated. 


tan  <*£ 


0 

0 


This  indeterminate  form  is  solved  by  applying  ^'Hospital's  Rule  to 

yield 


tan  aT 


“X1EX1E 

'*4E 


XlESin* 

cos  <f> 


-  X1E  *"an  ^  (50) 


or 


sin 


=  -  X1E  sin  <i> 


(51a) 


COS  Cfr 


cos  4> 


(51b) 


Since  sin  Ctp  and  cos  Op  are  just  the  negative  of  the  terms,  the  ap 
terminal  conditions  are  as  follows. 
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sin  Cfi 


=  Xlp  sin  0 


COS  Ct, 


=  cos  <p 


(52a 


(52b 


Linear  Development.  Once  again  the  terminal  surface  about 
the  pursuer,  P,  is  assumed  to  be  circular,  and  the  following  relation 
ship  represents  this  surface. 


where 


1 

2 


=  0 


Yj  =  R  sin  <p 


(53 


(54a 


and 

=  R  cos  <p  (54t 

with  the  angle  <p  representing  the  angle  around  the  terminal  surface 
measured  counterclockwise  from  the  pursuer's  local  horizontal. 

Taking  the  gradient  of  Eq.  (53)  yields  the  following  terminal 
relationships  for  the  linear  states  and  costates. 

V[Eq.  (53)]  = 
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The  remaining  conditions  are  determined  by  using  Eq.  (26)  for 
the  boundary  of  the  usable  part. 


H  =  v  •  X  =  AA  +  AX2Y2  +  AX^Y3  +  AX4y4 

ij  * 


f 


=  AX  jY^  +  AX2  2Yj  +  Y  j  +  Ug  sin  Otg  ~  Up  s 


iin  dpj 


»  ^  r 

AX3  -  Y2  +  UE  cos  Otj£  -  Up  cos  CCp  +  AX4  Y3  -  YJ  (56) 


This  equation  is  identical  to  Eq.  (39)*  From  Eq.  (25),  then,  Hjj(tf)  =  0, 
and  using  the  conditions  in  Eq.  (55),  this  relationship  reduces  to 

Y2  sin  <p  +  Y3  cos  <f>  -  R  sin  4>  cos  =  0 


or 


=  R  sin  <t>  -  Y2  tan  <t> 


(57) 


‘f 


To  evaluate  the  thrust  angle  (X  at  the  termina.  irface,  where 
Otj,  is  equal  to  Ctp,  it  can  be  seen  that 


tan  ctr 


AX2 


AX. 


0 

0 


By  applying  L'Hospital's  Rule,  the  value  of  a  at  the  terminal  surface 
is  determined. 


tan  a 


=  tan  <f> 


(58a) 
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(V 


aE  =  Ctp  =  0  at  tf 


(58b) 


Both  the  pursuer's  thrust  angle  and  the  angle  around  the  terminal 
surface  0  are  measured  from  the  pursuer's  local  horizontal.  When 
the  terminal  surface  has  been  reached,  then,  the  pursuer  would  desire 
to  have  maximum  velocity  outward  in  order  to  force  penetration  by  the 
evader,  thus  his  thrust  angle  would  equal  0,  The  evader,  still 
attempting  to  avoid  penetration,  would  obviously  point  his  thrust  in  the 
same  direction  as  the  pursuer's  thrust  to  neutralize  the  outward 
velocity  of  the  pursuer. 


!  C 
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IV .  Methods  of  Solving  the  Linear  Equations 

In  order  to  use  the  linear  equations  developed  in  the  foregoing 
chapters  'o  determine  the  barrier,  it  is  necessary  to  find  closed -form 
solutions  to  the  state  and  costate  equations.  This  chapter,  then, 
presents  the  methods  to  use  in  finding  these  solutions. 

The  Costate  Equations 

These  equations  are  to  be  used  in  backward  integration,  as 
mentioned  previously  in  Chapter  III.  In  order  to  do  this,  however, 
another  time  variable  must  be  introduced  which  corresponds  to  t£,  the 
final  time  at  the  terminal  surface.  Define  this  backwards  time  as 
T£  =  0  at  tf.  In  order  to  integrate  backwards,  the  signs  in  the  costate 
equations  are  reversed  so  that  proper  values  may  be  obtained  as  the 
trajectory  proceeds  outward  from  the  terminal  surface.  The  linear¬ 
ized  costate  equations  are  repeated  below  with  initial  conditions  at 
Tj  =  0  and  with  signs  reversed  for  backward  integration. 

Y  i 

AA1,=  AX2-AX4  AAj  (Tj=  0)  =  —  =  sin  0  (59a) 

AA2' =  AAj  -  AA3  AA2  (0)  =  0  (59b) 

AX3' AA3(0)=0  (59c) 

$ 

AAqt:  0  yielding  AA^  (Th)  =  a  constant  =  AA^(O)  =  cos  ^  (59d) 

The  solution  to  these  equations  is  approached  through  AA2  by  taking 
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another  derivative  with  respect  to  T  in  order  to  eliminate  terms  other 
than  AX2  and  constants.  This  yields  the  following  homogeneous 
solution. 

AX2h  =  A  sin  Tj  +  B  cos  T j  +  C  (60) 

By  differentiating  this  expression  twice  and  substituting  back  into  the 
differential  equation  for  AX2,  the  constant  C  may  be  found.  By  apply¬ 
ing  the  initial  conditions,  constants  A  and  B  can  also  be  found  yielding 
the  closed-form  equation  for  AX2  as  a  function  of  Tj. 

AX2  (T2)  =  sin  <f>  sin  Tj  +  2  cos  <f>  cos  Ti-  2  cos  <t>  (6la) 

By  applying  Eq.  (6la)  to  Eqs.  (59a)  and  (59c),  the  remainder  of  the 
closed-form  AX  equations  may  be  found. 

AXi  (Tj )  -  -  sin  <ft  cos  Tj  +  2  cos  <f>  sin  -  3  Tj  cos  (ft 

+  2  sin  <ft  (6 lb) 

AX3  (Tl)  =  -  2  sin  <f>  cos  Tj  +  4  cos  (ft  sin  Tj  -  3  cos  <f> 

+  2  sin  <f>  (6lc) 

AX4  (Tj)  =  constant  =  cos  <ft  (6ld) 

The  State  Equations 

The  state  equations  are  solved  in  the  same  manner  as  were  the 
costate  equations.  The  state  equations  with  initial  conditions  at  Tj  =  0 
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and  with  signs  reversed  for  backward  integration  are  repeated  below. 

In  these  equations  (f>  again  represents  the  angle  around  the  terminal 
surface  measured  counterclockwise  from  the  pursuer's  local  hori¬ 
zontal. 

y/ =-Y2  Yj(Tj  =  0)  =  R  sin  <f>  (62a) 

Y2  =  -  2Y3  -  Yj  -  (UE  -  Up)  sin  a  (62b) 

Y3'  =  Y2  -  (UE  -  Up)  cos  Of  Y3(0)  =  R  sin  $  -  Y2(0)  t  n  (62c) 
Y 4  =  -  Y3  +  Yx  Y4(0)  --  R  cos  <f>  (62d) 


The  closed-form  solution  to  these  linearized  state  equations  is 
approached  through  Y2;  however,  to  handle  these  equations,  approxi¬ 
mations  for  both  sin  Cl  and  cos  ct  are  required.  The  development  of 
these  approximations  is  presented  in  Appendix  A  with  the  following 
Taylor  Series  expressions  in  Tj  for  sin  Ct  and  cos  a  resulting. 


sin  Cl  =  + 


T1  + 


33C  _  3B3' 
2A3  8A5, 


(63) 


where 


A  =  sin  0 
B  =  2  sin  (f>  cos  <p 

C  =  -  —  +  2  cos^  + 
6 


si n^<p  cob^ 
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(64) 


where 

D  -  cos  £ 

E  =  2  sin  <t>  cos  0 

F  =  —  -  2  cos ctfa  -  sin  q>  cos  p 

The  approximations  for  sin  £1  and  cos  Cl  are  checked  for 
accuracy  in  the  equivalence  section  as  part  of  the  linear  to  non-linear 
comparison.  There,  when  the  negative  sign  is  applied  to  both  approx¬ 
imations,  the  resulting  trajectories  are  found  to  be  equivalent  out  to 
five  decimal  places  up  to  a  value  comparable  to  twenty  seconds  off  the 
terminal  surface.  In  Appendix  A  these  approximations  are  presented 
in  a  different  form  in  order  to  more  easily  determine  the  solution  to 
the  state  equations.  By  taking  another  derivative  of  with  respect 
to  T,  the  homogeneous  solution  is  found  to  be 


Y 


'H 


P  sin  Tj  +  Q  cos  Tj 


(65) 


Choosing  a  particular  solution, 

Y2p  =  K,  +  K2  Tj  +  K,  T,  2  +K4Tj3 


(66) 


The  total  solution  for  Y2  is  then 
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Y2  =  P  sin  Tj  +  QcosTj  +  Kj  +  +  K4Tj3  (67) 


Allowing  the  difference  in  thrust  constants  to  equal  (3  ((3  =  Ujr  -  Up) 
and  applying  the  initial  conditions,  the  constants  may  be  evaluated. 
Rather  than  list  these  constants  here,  however,  the  entire  solution 
is  included  in  Appendix  A. 

By  using  the  solution  for  Yz,  the  remainder  of  the  state  equa¬ 
tions  can  be  solved  in  closed-form  with  the  constants  again  evaluated 
by  applying  the  initial  conditions  at  the  terminal  surface.  The  closed- 
form  solutions  to  the  state  equations,  then,  are  as  follows. 


Yj  =  Pcos  Tj 


Q  sinT. 


1 


K1T1 


k3ti3 


K47y 


+  MM  (68) 


Y^  =  -  P  cos  Tj 


+  Q  sin  Tj  +  (Kj  +  (3A z)T1  +  (K2  +  (3BZ) 


+  (K3  +  (3cz) 


+  (K4  +  &DZ 


+ 


N 


(69) 


2  3 

Y4  =  2  P  sin  Tj  +  2Q  cos  Tj  -  (2K1  +  0AZ)  Zj_  -  (2K2  +  0B  z)l£ 


T  4  x  5 

(2K3  +  0CZ)  _1_  -  (2K4  +  /3DZ)  -i-  +  (MM-N)  Tj  +  S 


(VO) 


Once  again  the  constants  are  evaluated  in  Appendix  A  along  with  the 
entire  solution  to  these  equations. 
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V ■  Analysis  of  the  Barrier 

To  analyze  the  barrier,  trajectories  must  be  generated  out¬ 
ward  from  the  terminal  surface  to  provide  points  necessary  for 
constructing  the  barrier  itself.  The  method  used  in  this  thesis  to 
generate  these  trajectories  is  that  of  backward  integration.  This 
method  is  easy  to  apply  and  results  in  numerous  trajectories  from  the 
terminal  surface.  Although  the  analysis  of  these  many  trajectories 
may  be  tedious,  backward  integration  avoids  the  difficulties  of  iteration 
methods  while  providing  ample  information  on  the  barrier. 

In  applying  the  backward  integration  method,  two  computer 
programs  are  required,  one  for  the  linear  and  one  for  the  non-linear 
equations.  In  order  to  use  the  closed-form  solutions  to  the  linear 
equations  in  analyzing  the  barrier,  the  linear  trajectories  must  be 
equivalent  to  those  generated  by  the  non-linear  computer  program. 

[The  area  of  equivalence  will  be  discussed  in  a  later  section  of  this 
chapter.]  After  these  trajectories  are  compared,  linear  to  non-linear, 
the  points  on  either  side  of  the  barrier,  using  the  linear  trajectories, 
must  be  analyzed  to  determine,  where  possible,  whether  they  lie  in 
the  escape  zone  or  the  capture  zone.  This  chapter,  then,  presents 
the  linear  and  non-linear  computer  programs,  the  equivalence  analysis, 
and  the  methods  of  determining  whether  escape  or  capture  occurs  from 
any  given  starting  point. 
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The  Linear  Program 

In  developing  the  computer  program  to  nandle  tne  linear  state 
equations  in  closed  form,  a  primary  concern  was  the  data  to  be  used. 

In  Chapter  II  the  equations  of  motion  of  a  satellite  in  orbit  around  the 
earth  were  linearized  and  normalized  about  a  circular  reference  orbit. 
For  this  thesis  the  altitude  of  this  reference  orbit  above  the  earth's 
surface  was  selected  as  200  nautical  miles,  a  height  used  in  many 
manned  orbital  flights.  On  converting  to  the  earth's  canonical  units, 
the  radius  tq  of  this  orbit  is  1.  0581  distance  units  (DU)  with  a  refer¬ 
ence  velocity  Vq  of  0.  97216  DU/TU,  where  TU  is  the  earth's  canonical 
time  unit. 

Since  the  initially  specified  variable  Y£  is  the  difference 
between  the  pursuer's  and  evader's  radial  velocities,  which  have  been 
linearized  and  normalized  about  the  200  nautical  mile  reference  orbit, 
the  given  value  of  Y^iO)  at  the  terminal  surface  will  have  to  be  quite 
small.  In  ordex  to  gain  equivalence  between  the  linear  and  non-linear 
programs,  an  arbitrary  value  of  Y2(0)  was  taken  as  0.  005.  This  value 
is  equivalent  in  size  to  a  relative  velocity  difference  of  approximately 
130  feet  per  second.  The  value  of  Y2(0),  however,  could  actually  be 
varied  while  holding  the  other  variables  constant  as  an  alternative  in 
seeking  more  information  about  the  shape  of  the  barrier. 

In  conjunction  with  the  arbitrary  selection  of  Y2(0)  is  the  selec¬ 
tion  of  the  value  for  the  difference  in  normalized  thrusts,  Uj-  -  Up, 
labeled  BT  in  the  computer  program.  The  numerical  value  of  BT  was 
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selected  to  give  the  pursuer  a  slight  advantage  since  equal  values  of 
UE  and  Up  would  make  DT  =  0,  thus  removing  the  control  from  the 
solution.  BT,  then,  was  given  the  value  -  0.  04  reflecting  the  fact  that 
the  pursuer's  thrust  is  greater  than  the  evader's  thrust.  This  value 
of  normalized  thrust  difference  resulted  from  applying  thrusts  of  a 
magnitude  equivalent  to  providing  just  over  one-tenth-g  of  acceleration. 
With  the  value  of  BT  specified,  only  the  radius  of  the  terminal  sur¬ 
face  R  and  the  angle  <f>  on  that  surface  remain  unspecified. 

Since  the  value  of  <f>  will  obviously  vary  around  the  terminal 
surface  in  establishing  the  boundary  of  the  usable  part,  this  term  will 
remain  as  a  variable  to  be  changed  readily  in  both  programs.  Since 
the  radius  of  the  terminal  surface  R  is  a  non -dimensional  quantity, 
its  value  must  be?  jelected  in  reference  to  the  siz.  the  reference 

orbit  radius.  <  ike  radius  of  the  terminal  surface  ere  equivalent  to 
1000  ft,  the  value  of  the  non-dimensional  R  would  be  4.  78  x  10"J. 

With  these  vc-ues  of  <f>  and  R,  the  initial  conditions  for  the  linear  pro¬ 
gram  were  completely  specified.  These  initial  conditions  were  then 
applied  to  the  closed-form  solutions  to  the  linear  equations  presented 
in  Chapter  IV  to  generate  the  required  trajectories  backward  in  time 
from  the  terminal  surface.  A  sample  of  this  computer  program  is 
included  in  Appendix  B. 

The  Non-linear  Program 

For  equivalence  between  the  linear  and  non-linear  programs. 
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data  values  comparable  to  those  in  the  linear  program  must  be 
selected.  Here,  however,  eight  initial  values  are  necessary  in  order 
to  completely  specify  the  initial  conditions  due  to  the  dimensionality  of 
the  non-linear  program,  whereas  only  four  initial  conditions  were 
required  in  the  linear  program.  The  values  for  Tq  and  Vq  are  not  con¬ 
sidered  as  initial  values  since  these  are  only  used  as  references  and 
are  not  specifically  used  in  the  programs.  In  determining  these  initial 
values,  first  of  all  consider  the  values  that  are  common  to  both  pro¬ 
grams. 

The  values  of  R,  <f>,  and  BT  would  be  the  same  as  in  the  linear 
program  with  BT  split  into  individual  values  of  =  0.  1 1  and  Up  =  0.  15 
giving  the  same  difference,  BT.  The  other  common  value  is  that  of 
Y£(0),  here  representing  again  the  difference  between  the  players' 
radial  velocities,  assuming  that  the  non-linear  relative  difference  is 
approximately  equal  to  the  linear  relative  difference.  The  other  values 
for  this  program  are  selected  mostly  from  considering  the  previously 
given  values. 

The  remaining  values  at  termination  to  be  specified  are  the 
pursuer's  state  variable  initial  conditions.  Since  the  evader's  initial 
conditions  may  be  determined  from  these  values,  the  pursuer's  values 
are  the  only  ones  required.  Eqs.  (47)  and  (49)  for  the  relationship  of 
these  initial  conditions  are  repeated  below  for  easy  reference. 
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X1E  -  Xlp  =  R  sin  (47) 

X4E  '  X4P  =  R  cos  * 

(X2E  -  X2P>  sin  *  +hr^-  •  cos  ^  =  0  (49) 

\X1E  X1P^ 

Since  these  values  are  to  be  determined  at  termination,  the  pursuer's 
values  are  arbitrarily  assumed  to  approach  the  values  of  the  reference 
orbit  for  easier  handling.  Therefore,  the  values  of  the  normalized 
radius  and  tangential  velocity  will  be  one,  and  the  value  of  the  normal¬ 
ized  radial  velocity  will  be  zero.  This  leaves  the  angle  as  the  only 
unspecified  value. 

In  selecting  dp,  the  value  of  =  A#E  -  from  the  linear 
program  was  considered  since  it  is  presumed  that  Y^(0)  will  be  equiva¬ 
lent  to  the  difference  between  0E  and  $ p  from  the  non-linear  program. 
Referring  to  Eqs.  (47)  and  remembering  that  the  magnitude  of  R  is 
4.  78  x  10“^,  the  value  of  R  cos  <p  and  hence  the  difference  in  6^  and 
dp  will  be  quite  small.  The  value  of  #p  is  then  arbitrarily  chosen  to 
be  4  x  10  .  This,  then,  completely  specifies  the  data  values  for  the 

non-linear  program.  Once  again  a  sample  of  this  program  is  pre¬ 
sented  in  Appendix  B  for  reference  or  future  work. 

Equivalence 

Once  the  data  values  for  both  computer  programs  were  estab¬ 
lished  as  of  comparable  magnitude,  the  trajectories  produced  by  each 
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program  were  compared  for  equivalence.  The  search  for  equivalence 

started  with  attaining  the  proper  values  specified  previously  and  with 

finding  the  proper  step  size  for  the  non-dimensional  time,^j.  This 

quantity's  numerical  size  was  determined  from  the  relationship, 

V  t 

= _ 2_  >  where  t  is  in  time  units.  Since  the  canonical  time  unit  is 

ro 

so  large,  the  value  of  Tj  was  successively  reduced  by  a  factor  of  one- 
tenth  in  searching  for  equivalence.  At  a  step  size  of  0.  001  the  tra¬ 
jectory  values  from  both  the  linear  and  non-linear  programs  came  to 
within  2%  of  each  other  out  to  an  equivalent  of  two  hundred  seconds  off 
the  terminal  surface.  Some  of  the  values  that  were  compared  for  the 
angle  <f>  of  355  degrees  are  listed  in  Table  I. 

Both  the  linear  programs,  one  for  integration  of  the  linear  state 
and  costate  equations  and  one  for  the  closed-form  solution  equations, 
were  found  to  be  equivalent  out  to  a  comparable  value  of  at  least  two 
hundred  seconds  off  the  terminal  surface.  The  importance  of  gaining 
equivalence  for  the  closed-form  equation  trajectories  lay  in  the  desire 
to  use  these  closed-form  equations  for  analyzing  the  starting  point 
locations.  This  analysis  provides  the  final  section  in  this  chapter  on 
analyzing  the  barrier. 

Starting  Point  Locations 

The  analvsis  of  the  harrier  was  ronrluded  hv  a  searrh  for  a 

J  ‘  /  ' 

method  of  determining  whether  a  given  starting  point  lay  in  the  escape 
zone  or  in  the  capture  zone.  This  problem  had  to  be  approached 
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Table  I 

Equivalence  Values,  <f>  =  355  Degrees 


Non-linear  Integration  Values 


Y(l) 

V  (2) 

Y(3) 

Y  (4) 

TAUj 

-4.  16629E-06 

5.  OOOOOE-03 

4.  33303E-04 

4.  7b  18 IE-05 

0. 

-9.  16400E-06 

4.  99549E-03 

4.  78149E-04 

4.  71557E-05 

0,  001 

-5.  39272E-05 

4.  95020E-03 

8.  S1359E-04 

4.  07532E-05 

0.  01 

-1.  27482E-04 

4.  85243E-03 

1.  55153E-03 

2.  1 1428E-05 

0.  025 

-2.  46 146E-04 

4.  62799E-03 

2.  66208E-03 

-3. 62128E-05 

0.  05 

-3.  58246E-04 

4.  32742E-03 

3.  76253E-03 

- 1 . 24075E-05 

0.  075 

-4.  61890E-04 

3.  95161E-03 

4.  85054E-03 

-2.  41966E-05 

0.  1 

-5.  55208E-04 

3.  50148E-03 

5.  92382E-03 

-3.  89311E-04 

0.  125 

-6.  36353E-04 

2.  97796E-03 

6.  98012E-03 

-5.  65437E-04 

0.  15 

-7.  03502E-04 

2.  38203E-03 

8.  01 727E-03 

-7.  69576E-04 

0.  175 

-7.  54898E-04 

1. 71466E-03 

'9.  03313E-03 

-1.  00087E-03 

0.  2 

Linear  Closed-Form  Equation  Values 


Y  (1 ) 

Y  (2) 

~  Y(3 ) 

Y  (4) 

TAUj 

-4. 16629E-06 

5.  OOOOOE-03 

4. 

33303E-  04 

4.  76  181E-05 

0. 

-9. 16410E-06 

4.  99559E-03 

4. 

78147E-04 

4.  7 1 557E-05 

0.  001 

-5. 39281E-05 

4.  95027E-03 

8. 

81361E-04 

4.  07532E-05 

0,  01 

-1.  27482E-04 

4.  85223E-03 

1. 

55160E-03 

2.  1 1405E-05 

0.  025 

-2.46126E-04 

4.  62626E-03 

n 

u. 

66231E-03 

-3.  62340E-05 

0.  05 

-3.  58 148E- 04 

4.  32272E-03 

3. 

76275E-03 

-1.  24141E-05 

0.  075 

-4.  61630E- 04 

3.  94338E-03 

4. 

85025E-03 

-2.  42  1 00E-05 

0.  1 

-5.  54709E-04 

3.49109E-03 

5. 

9222 1 E- 03 

-3.  89519E-04 

0-  125 

-6.  35609E-04 

2.  96983E-03 

6. 

97608E-03 

-5.  65694E-04 

0.  15 

-7.  02668E-04 

2.  38468E-03 

8. 

00946E-03 

-7.  69817E-04 

0.  175 

-7.  54362E-04 

1.  74182E-03 

9. 

02005E-03 

-1.  00098E-03 

0.  2 
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through  the  linear,  closed-form  solution  equations  as  mentioned 
previously.  Since  the  barrier  is  defined  in  terms  of  four  dimensions 
in  three  parameters,  the  method  of  specifying  starting  point  locations 
had  to  be  completely  analytical  through  use  of  the  linear  equations. 

Also,  since  time  limitations  precluded  any  extensive  examination  of 
this  method,  this  development  really  results  in  a  theoretical  method 
for  considering  given  starting  points  near  the  barrier. 

The  approach  to  discovering  a  method  to  locate  given  staicing 
points  in  relation  to  the  barrier  was  to  find  a  normal  to  the  barrier 
surface  which  would  pass  through  the  given  point.  In  three-dimensional 
surfaces  the  normal  is  based  on  the  orthogonality  of  the  coordinate  sys¬ 
tem  in  which  the  surface  is  presented.  Since  the  four  coordinates  of 
this  thesis  are  also  orthogonal,  the  three-dimensional  normal  theory 
was  extended  to  handle  the  four -dimensional  barrier  surface  that  was 
found  in  this  study. 

The  normal  to  a  four-dimensional  surface  was  formulated  from 
the  three-dimensional  development  of  R.  C.  Buck  (Ref  6:335).  In  this 
reference  the  normal  to  a  three-dimensional  surface  as  a  function  of 
two  parameters,  u  and  v,  was  found  as  follows. 

t 

To  extend  this  equation  from  three  dimensions  in  two  parameters  to 
four  dimensions  with  three  parameters:  <p ,  Y^(0),  and  Tj,  the  follow¬ 
ing  normal  equation  was  proposed. 
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N  = 


d(Y2,  Y3>Y4)  >  d(Y3 ,  Y4,  Y i ) 
d(<?»,Y2{ 0),rj)  '  d(^Y2(0),Ti) 


d(Y4,  Yi,  Y2) 
<*(0.  Y2(0),t!) 


rf-.ri-  y2)  y3) 
d’  ■'  ^2<°)»rl) 


(71) 


In  this  equation  the  Y-terms  refer  to  Eqs.  (67)  through  (70),  the  linear, 
closed-form  solution  equations.  By  evaluating  the  four,  three-by- 
three  determinants  in  Eq.  (71),  the  normal  at  any  point  on  the  barrier 
may  be  determined.  Due  to  the  complexity  of  these  determinants,  it 
was  not  deemed  beneficial  to  list  the  final  equation  form  for  this  four¬ 
dimensional  barrier.  The  easiest  method  of  dealing  with  the  normal 
would  be  to  take  the  partial  derivatives  of  the  Y-equations  with  respect 
to  the  three  parameters  and  evaluate  these  at  a  point  along  the  barrier. 
Once  all  the  partial  derivatives  for  a  given  point  have  been  evaluated, 
the  determinants,  and  thus  the  normal  itself,  could  be  determined. 

Once  .he  normal  has  V  en  determined,  the  obvious  question 
deals  with  wnerc  the  given  point  o' 1  the  barrier  i  a  relation  to  the 
barrier-  On  generating  normals  all  along  the  barrier,  one  will  surely 
pass  through  the  giver,  starting  point  or  point  exactly  opposite  to  the 
direction  through  the  given  point.  By  taking  the  dot  pr  oduct  between 
any  normal  vector  and  the  vector  from  the  base  of  that  normal  to  the 
given  point,  the  angle  between  the  two  vectors  may  be  found.  By 
searching  in  the  four  coordinate  directions,  it  should  be  possible  to 
find  the  point  on  the  barrier  where  the  normal  passes  through  the  given 
starting  poi.it  or  180  degrees  away  from  the  point.  Knowing  this  would 
tell  whether  the  point  was  on  the  positive  or  negative  side  of  the  barrier 
and,  therefore,  in  the  capture  or  escape  zone  respectively. 
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VI.  Conclusions  and  Recommendations 

Conclusions 

The  objective  of  this  study  was  to  determine  a  differential  game 
barrier  for  a  pur  suit -evasion  game  between  two  spacecraft  in  orbit. 
Throughout  the  chapters  of  this  thesis,  this  general  area  of  concern 
was  broken  down  into  smaller  sections  dealing  with  the  search  for  this 
barrier.  These  sections  dealt  with  the  development  of  the  linear  and 
non-linear  state  and  costate  equations,  the  theory  of  differential  game 
barriers,  and  finally  an  analysis  of  the  barrier  itself. 

The  state  and  costate  equations  were  used  to  generate  trajec¬ 
tories  backward  in  time  from  the  terminal  surface.  These  trajectories 
were  then  compared,  linear  to  non-linear,  and  once  equivalence  had 
been  verified,  the  barrier  for  this  problem  was  determined.  However, 
this  barrier  was  a  four-dimensional  surface  that  depended  on  three 
parameters:  <f>,  YgCO),  and  r^.  With  =  o  at  the  terminal  surface, 
one  specified  set  of  values  for  <f>  and  Y2(0)  established  a  unique  bound¬ 
ary  of  the  usable  part  in  four  space.  This  is  obviously  difficult  to 
demonstrate  graphically;  however,  by  making  a  ratio  of  the  relative 
velocities,  Y2/Y3,  a  three-dimensional  presentation  of  the  barrier  is 
possible.  Figures  4  through  7,  then,  portray  the  boundary  of  the 
usable  part  in  the  four  quadrants  of  a  cylindrical  terminal  surface. 
Along  with  the  boundary  of  the  usable  part  in  Fig.  4  is  a  portion  of  the 
barrier  itself  out  to  a  value  equivalent,  to  one  and  one-half  seconds  off 
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the  terminal  surface.  To  better  present  this  barrier,  two  planes  were 
passed  through  the  cylinder  parallel  to  the  plane  containing  Yj  and  Y4 
at  values  of  Y2/Y3  equal  to  plus  and  minus  one  hundred  in  the  scale  of 
Fig.  4.  These  two-dimensional  barrier  segments  out  to  an  equivalent 
of  twenty  seconds  off  the  terminal  surface  are  shown  in  Figures  8 
and  9. 

Once  the  barrier  was  determined,  then,  meeting  the  main 
objective  of  this  study,  the  barrier  itself  was  analyzed  to  determine 
the  location  of  various  starting  points  lying  just  off  the  barrier.  A 
theoretical  presentation  of  a  method  to  determine  these  locations  was 
presented  in  the  last  chapter.  Time  limitations,  however,  precluded 
a  complete  evaluation  of  this  method,  as  mentioned  previously;  there¬ 
fore,  this  portion  of  the  objective  was  only  partially  completed. 

Recommendations 

The  first  and  most  obvious  recommendation  would  be  to  extend 
the  analysis  of  the  barrier  to  determine  the  location  of  starting  points 
in  the  playing  space  lying  just  off  the  barrier.  Since  the  linear, 
closed-form  equations  have  been  checked  for  accuracy  with  the  non¬ 
linear  equations,  these  linear  equations  provide  a  ready  starting  point 
for  this  further  analysis.  Once  the  technique  for  determining  the 
normal  to  a  four-dimensional  barrier  surface  has  been  evaluated,  the 
location  of  the  given  starting  points  could  be  determined. 

Another  area  for  further  study  would  be  to  extend  the  analysis 
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Barrier  Cut  at  Y7/ Y,  =  +  1.  00 


Mw  'uA/  f jA-j 

of  the  barrier  in  three  dimensions.  By  making  more  cuts  along  the 
Y2/Y3  axis,  the  barrier  could  be  completely  defined  allowing  graphi¬ 
cal  analysis  along  with  the  jbove  mentioned  analytic  analysis.  In 
addition  to  this  recommendation  would  be  one  to  vary  the  value  of  Y2(0)» 
the  relative  velocity  difference  in  the  radial  direction  at  the  terminal 
surface.  This  would  generate  a  changing  barrier  surface  which  could 
provide  information  concerning  the  velocity  requirements  fo*  either 
capture  or  escape. 
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Appendix  A 


Sine  and  Cosine  Approximations 

In  order  to  use  the  state  equations  effectively,  it  is  necessary 
to  get  the  sin  Ct  and  cos  d  terms  into  expressions  that  can  be  solved. 
This  requires  expanding  the  expressions  for  sin  d  and  cos  Ct  in  Taylor 
Series  expansions  about  the  point  =  0,  where  T ^  is  the  non- 
dimensional,  backward  integrating  time.  In  the  linear  equations,  as 
stated  in  Chapter  III,  the  optimal  control  angles  for  the  pursuer  and  the 
evader  were  found  to  be  equal;  therefore,  only  two  approximations,  one 
for  the  sine  and  one  for  the  cosine,  are  necessary.  Since 


sin  Ci 


A  A  2 


V‘ 


Aa|  +  AA: 


(72) 


it  was  necessary  to  first  square  this  expression  before  expanding.  This 
results  in  an  equation  of  the  form 


sin^  Cl 


**2 

A  A  2  +  A  A  ^ 


*  (0»T!) 
g  (<£>Tl) 


(73) 


This  is  readily  expanded  in  the  following  manner,  assuming  that  is 
small. 
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~dl 

IJJL 

‘dzi 

2  dl  dz 

sm2a  =  f cp.ro 

4* 

dT  1 

.  dT\ 

T,  +  I 

dT? 

dTx  dT\ 

g(^*Tl) 

O 

11 

.8 

71 

n-»  2 

_  g 

g2 

i  JjL.  2f (dj^ 


dT, 


dr , 


3  ‘ 


(74) 


Ti  =0 


J  '1 


where  ^ 

m>rx)  =  s s  in^Tji  +  4  cos^^cos^T^  +  4cos2^ 

+  4  sin  <f>cos  <p  sinT^  ccsT^  -  4  sin  cos  0  sinTj  -  8  cos^cos^  (75a ) 

and 

g^Tj)  =  4  +  4  cos^  +  12  sin  cos<j6  sir/T^  -  12  sin^>  cos  ^>sin‘7jcos‘7| 

+  sin^sin^  -  8  cosTJ  +  l6cos^sin^  t  97J2cos2$ 

+  127^  sin<p  cos^  cosTj  -  24Tj  cos20  sinTj  -  1 2*7^ sin cos <f>  (75b) 


At  Tj  =  0,  these  terms  reduce  to  zero..  By  taking  derivatives  of  f  and  g 
with  respect  to  7^  and  then  evaluating  at  7^  =  0,  the  following  expres¬ 
sions  result. 


di_ 

dTx 


=  0 


Tr° 


dzS 


=  2  sin^0 


Tr° 


(75c) 


(75d) 
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lL 

<*Ti' 


-  -12  sin  <f> cos 


lrx=0 


d4f 
d  rj4 


-8  +  32  cos*  (ft 


rl=° 


48- 

dr i 


=  0 


!r  j  --  0 


(pg 


dr ^ 


=  2 


r1=0 


O' 


d3g 


dr j3 


d4g 

^7 


=  0 


r,=0 


16  -  24  cos  <f> 


rl=0 


(75e) 


(75f) 


(75  g) 


(75h) 


( 7  5 ;  > 


(75j) 


The  first  term  in  Eq.  (74)  is  then 


g(fcr!> 


0 

0 


Applying  L'Hospical's  Rule  twice  to  this  indeterminate  expression 
yields  the  following  relationship. 


d2f 


4?S- 

d  2 


sin^ 


ri=0 


(76) 
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This  then  is  the  first  coefficient  of  Eq.(74).  The  second  coefficient  is 
approached  in  the  same  manner  after  putting  the  expression  over  a 
common  denominator. 


g  d£_  -  f  dg  ' 
dr j  d  rj 

~2 


(77) 


Applying  L'Hospital's  Rule  four  times  finally  results  in 


and 


on 


g  di  _  f  dg 


dr. 


dri 


_  2  d3f 
drl2  dr1^ 


drl"* 


=  6 


.2 
±JL. 
dry 


(78a) 


(78b) 


At  r^  0,  Eq.  (78a)  equals  -  48  sin<£  cos  0  and  Eq.  (78b)  equals  24  to 
yield  the  second  coefficient  for  Eq.  (74). 


df 

dr  i 


-f  -Jf 
d  1 


g 


2  sin  </>  c.os  ^ 


-*  r1=0 


(79) 


The  third  coefficient  of  Eq.  (74)  is  also  put  over 


a  common  denominator 


giving  this  expression. 


! 

~  ,rf2  f 

wf 

r/  <r 

*2,. 

1  | 

|g  drl* 

"  2g  — - 

drt 

ef  +  2f 

2 

L 

(80) 
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By  applying  L'Hospital's  Rule  six  times  as  before  and  evaluating  at 
=  0,  the  third  coefficient  of  Eq.  (74)  is  determined  to  be 


1 

2 


dZl 


2  dtdg  i 

d^d  Tj 
— 


2 

d  g 
dT1 2 

“F" 


2f(Vg 

<nri 


g' 


7^0 


5.  +  2  cos ^d>  +  sin^<£cos2<6 
6 


(81) 


Substituting  the  results  of  Eqs.  (76),  (79),  and  (81)  into  Eq.  (74)  yields 
the  following  expression  for  sin  Ct. 


sin  Ct  =  +  ^sin^0  -  (2  sin^cos$)7j  +  -g-  +  2cos^fsin^cos^|Tl^j 


1/2 


(82) 


Since  this  expression  is  actually  the  square  root  of  the  function, 
another  approximation  is  required.  The  above  expression,  then,  is 
expanded  in  a  Taylor  Series  expansion  about  Tj  =  0  as  follows,  the 
coefficients  of  Tj  being  lumped  into  constants  for  ease  of  handling. 


F(T,)  =  (A*  -  BT  +  CT,  ) 


2.1/2 


F(0)  =  A 


F'(Ti)  =  7  (A2  -  B7i  +  CT12)”1/2  (-B  +  2CTX)  F'  (0)  =  ^ 

«  ZA 
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f"  ( r j)  =  -  -  (A2  -  Btj  +  CTi2)"3/2  (-B  +  2CTJ)2  4  C  (A2  -  BTj 


4  CT,  ) 


2'-1/2  f"  (0)  =  ,  C 


4AJ  A 


F"'  (Tx)  =  |  (A2  -  B Ti  4  CTl2)"5^2  (-B  4  acTj)3  -  (A2  -  BTj 

4  CTi2)_3/2  (-B  4  ZcTjMC) 


-  §  (A2  -  BTx  4  CTl2r3/2  (-  B  4  2CTJ) 


F"'  (0)  = 


3BC  _  3B- 
2AJ  8A~ 


Therefore, 


|  f,  B  T  ,1c  B2  \T12  ,  f 3UC  3B3\Ti3l 

«  -  i  [A  '  Iati  p  -  4Xr]  nr  3  -  il?]— J 


where 


.A  =  sin  (f> 

B  =  2  sin  </>  ccs  </> 


C  =  -  —  4  2  cos2^>  4  sin2<jScos20 


(83a) 


(83b) 

(83c) 


(83d) 


In  a  similar  manner  as  abcve,  the  cos  C(  approximation  is 
determined;  however,  in  this  case  the  identity  cos  Ct  =  1  -  sin  Ct  is 
used.  From  this,  an  expansion  of  the  square  root  is  accomplished  as 
before  wielding  the  following  cos  Ct  annroximation- 


S  a  =  +  [d  +  I,  Tj  +fe  -  tfev-  - 

"  p  2D  |D  4D3/  2  18DJ  2d?J  o  J 


(84a) 
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where 


D  =  cos  <f> 

(84b) 

E  =  2  sin  <f>  cos  <f> 

(84c) 

F  =  —  -  2  cos^  d>  -  sin^  <t>  cos^d> 

6 

(84d) 

Linear  State  Equation  Solution 

The  linear  state  equations  from  Chapter  III  are  repeated  below 
for  easy  reference  with  the  signs  reversed  for  backward  integration 
and  with  the  determined  initial  conditions. 


C4 

>-< 

1 

It 

r-H 

>< 

Yj^  =  0)  =  R  sin<£ 

(85a) 

Y2'  =  -  2Y 3  -  Yx  -painCC 

(85b) 

Y,'=  y2  - /3cos  a  y3(0) 

=  R  sin  <f>  -  Y2(0)  tan  <f> 

(85c) 

Y4,=  '  Y3  + 

Y^(3)  =  R  cos  0 

(85d) 

In  these  equations,  P  =  Up  -  Up,  the  difference  in  normalized  thrusts, 
is  the  non-dimensional  time  for  backward  integration,  and  fp  is  the 
angle  around  the  terminal  surface  measured  counterclockwise  from 
the  pursuer's  local  horizontal.  Rather  than  immediately  use  the 
approximations  for  sin  Ot  and,  cos  Ct  from  the  first  section  of  this 
appendix,  a  new  set  of  constants  is  selected.  In  making  this  selection 
the  proper  sign  of  Eqs.  (83a)  and  (84a)  had  to  be  chosen.  Only  after  the 
equivalence  analysis  was  completed  were  the  proper  signs  known; 
therefore,  the  following  expressions  are  negative  in  order  for 
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equivalence  to  occur  between  the  non-linear  and  linear,  closed-form 
computer  programs. 


where 


sin  a  =  -  (Aj  +  +  CiTi2  +  DjTi3) 

(86a) 

Aj  =  A 

(86b) 

Bj  =~ 

1  2A 

(86c) 

c  =  I  l£  -  sLj 

1  2(A  4A3| 

(86d) 

„  _  1  f  3BC  3B3\ 

1  6  IZA3  8A5) 

(86e) 

and 

cos  a  =  -  (Az  4-  B2Tj  +  CzT^  +  D2Ti3) 


(87a) 


where 


a2  = 

D 

(87b) 

b2  = 

E 

2D 

(87c) 

C2=ji 

F  E2  | 

D 

4D  • 

(87d) 

D24  1 

f  3E3  3EF  \ 
i8Db  2D3  / 

(87e) 

The  constants  A  through  F  come  from  the  first  section  of  this  appendix 
where  they  are  evaluated. 
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By  differentiating  Y2*  again  with  respect  to  T  and  substituting 
in  from  the  other  state  equations,  an  equation  relating  Y2  terms  and 
constants  results. 

Y2"+  Y2  =  -2p(Az  +  B2TJ  +  C2Ti2  +  D27i3)  +/3  (Bx  +  20^1 

+  SDfi2)  (88) 

or 

Y2  =  P  sinTj  +  Q  cosTi  (89) 

H 

Choosing  now  a  particular  solution  to  the  above  equation, 

Y2p  =  Kj  +  K2Tx  +  K3T^  +  K4T!3  (90) 

and  differentiating  with  respect  to  Tj  the  constants  through  K4  may 
be  evaluated. 

2K3  +  6K4Tj  +  Kj  +  K2Tl  +  K3Tx2  +  K4T j3  =  -/3[  (2A2  -  Bj) 

+  f2B2  -  +  (2C2  -  SDjjTi2  +  2D2Tj3]  (91) 


Kj  =  -  /3(2A2  +  6D1  -  Bj  -  4C2)  (92a) 

K2  =  -  iS(2B2  -  2CX  -  12D?)  (92b) 

k_  =  -  R12C-,  -  3D .  t 

'  £.  1' 

K4  =  -  2  (3d2  (92d) 
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This  yields  the  following  total  solution  for  Y2* 

2  3 

Y2  =  Y2  +  Y2p  =  P  sinT^  +  Q  cosT1  +  Kj  +  +  K 4TX  (93) 

By  use  of  the  initial  conditions  in  Eqs.  (85)  the  constants  P  and  Q  may 
be  determined. 

Y2  (Tj  =  0)  =  Q  +  Ky  (94) 

Therefore, 

Q  =  Y2(0)  +  /3  (2A2  +  6DX  -  B1  -  4C2)  (95) 

and 

Y2,(0)  =  P  +  =  -  3R  sin<0  +  2Y2(0)  tan  <p  -  sin<p  (96) 

or 

P  =  (-3R  -/S)  sin  0  +  2Y2(0)  tan0  +  )3(2B2  -  2Cj  -  12D2)  (97) 

where  Y2(0)  is  a  specified  initial  condition  for  this  thesis. 

By  using  Eq.  (93),  the  remainder  of  the  closed-form  solutions 
to  the  linear  state  equations  may  be  found, 

Y  |'  =  -  P  sin*T^  -  QcosT^  -  -  K2Tj  -  K37i2  -  K^Tj 3  (98) 

_  k7t2  k?t3  k/7*,4 

Y i  =  P  cos  Tj  -  Q  sinTj  -  KjTj  -  _i_L.  -  i,  -  _  4  i._  +  MM  (99) 

2  3  4 

where  MM  is  determined  by  use  of  the  initial  conditions  tc  be 

MM  -  (4R  +0)  sin  0  -  2Y2  (0)  tan  0  -  /3  (2B2  -  2Cj  -  12D2)  (100) 
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K7T2  K,T3  K47T4 
Y3  =  -  P  cos  Tj  +  Q  sinTj^  +  KjTj  +  -  ■■■  +  ■  .—  l- 


P[A2rI  +  5fl  +  ^£i!+^V) 


+  N 


where 


N  =  (  -2R  -#)  sin  <£  +  Y2(0)  tan  <?>  +  /3  (2B2  -  2CX  -  12D2) 

T  2  'T.  2 

Y4  =  2P  sinT!  +  2Q  cosTj  -  (2KX  +  /3A2)  _L_  -  (2K2  +  0B2)  — 


f  4  7*i  5 

-  (2K,  +  /3C?)  li_  -  (2K4  +  0D2)  -L-  +  (MM  -  N)  T,  +  S 
3  c  12  20  1 


where 


S  =  R  cos  <f>  -  2Y2(0)  -  2/3<2A2  +  6Dj  -  Bj  -  4C2) 


(101) 


(102) 


(103) 


(104) 
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Appendix  _B 


Non-linear  Integration  Program 

PROGRAM  NCNLIMIKPUT  ,CU T F IT  ,  T A FE5* I NPUT ,  T A F E6 =OUT P UT> 

Cl  PENSION  XE  (4)  ,XP(4),CX£(4)  ,CXP(4),LDA£(4>  5LQAP<4>*0L0AE<3)| 

1  CL  DAP  {  2)  ,xrAVG(4),L£AVC<2),XFAVG(4),LPAVG(2)  >  Y  <  4)  ,SXE  (  4)  , 

2  S XD ( 4 )  »S  EXE ( 4 ) tSOXF  (4) »  S  C  L  C  E  (3)  >SD  LOP  (  2  )  »SLE  (3) »  S IP ( 3 ) 

SEAL  LOAE  ,LCAPf LEAVG *LFAWC 

P£A0<5, 13 )  FEE,R,XP  d)  ,XF  <2>  ,>P(3)  ,  XP<4) 

13  FORMAT ( Ecll. 5) 

FEA3C5,11)  V",Pd,Y23,TAL,LE,lF 

11  FORMAT { fell. ?) 

READ  (5 , 12)  LOAE(2),LCA£(2),LCAF(2),LDAP(3i 

12  F  CCMAT { hE 11 • 5 ) 

WRxTE  (6.EP0) 

513  FORMAT  (MM 

13  F  r  3  =  (  FEE*3. 14159) /18Q. 

TANPMI  =  SIN  (FHI) /C  CS (PHI) 

LCAE(l)  =  SIN(PHT) 

L C Ac ( 4 )  =  COS  (PHI ) 

LCAP(l)  =  -SIM  PHI) 

LCAP14)  =  -COS ( PHI) 

WRITE  ( 6  »2C) 

2  J  FORMAT ( *  PHI  R  XP(1)  XP(2) 

1  y P ( 2 )  >F ( 4)  TAl  *) 

WRITE  IE, 21)  FEE,R,  XFU)  ,XF  (2)  ,XP(3)  ,XP  (4)  ,TAU 
21  FORMAT  (1P7E13. 5) 

CT  =  'l.  JuCOl 
K=C 
K  =  1 

X  E (1 >  =  R*SIN<PHI) +XF (1) 

XE  (2)  =  Y2',+XP  (2) 

XE  (3)  =  R*SIN(PH)~(XE  <  2  )  —  X  F  (2)  )*TANPHI  +  XF  (2) 

X  E (4 )  =  ovCOS (PHI )  +  XF (4) 

Y  (1) =XE  (l)-XP  (1) 

Y  (2)  =X  E  (2  )-X°  <2  ) 

Y  <  2 )  =  X  E  (2  )-XP  (3) 

Y  ( 4)  =  ( X  E  (  4 )  —XF  (  4)  ) 

FTS=SQDT (Y(l)¥¥2+Y (  4  )  * *2 ) 

VTCT=SQRT (Y(2)**2+Y  \2)**2) 

VRAT-Y ( 2) /Y ( 3) 

WRITE  <  6, 11C ) 

110  FORMAT  <*  Yd)  Y(?)  Y(3)  Y(4) 

1  Y  ( 2  !  / Y  ( 7)  VTCT  TAL'“5 

WRITE  <c,25)  Y  ( 1)  ,Y  (2) ,Y  (2 )  ,  Y  (4)  ,RTS,VRAT  ,VtCT,TAU 
25  FORMAT  ( 1P8E15,  5) 

1 n  IF  (TAU.GT.e.)  GC  TC  111 

CXE(2)=-(XE(3)**2/XE(1)-1./XE  (1 )  **  2-iJF*  X  £  d  )*SIN  CPU)  ) 
r>F(3)  =  -C-X«*(?t  *XF<  2  )/xr  d  MlE*COS<°HI)  ) 

CXF(?)=-(Xn(3),l<l2/XF(l)-l-/>Fd)  **2+UP*XR  d  )  M+SIN  (PHI  ))) 
CXFt3)=-(-XPC2»*XP(  2)/XF(J)4LFM+COS(PHI))) 

C-C  TO  112 
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1U  CXE(2)=-(XE  (  3 ) ** 2/X  E  Cl )  -  1  *  /  >  E  (i)  +  UE*  L  C  AF  (2)/($QPT(L0AE(2)**2 

1  +LOAE (  3)  **2)  ) ) 

-  CXEC3)  =- t-XE  (2)  *XE(  3)/XE  (1  >  m^LOAE  (?>/ (SCR  7 (LDAP ( 2)**2+LD  A£ ( 3) 
1  **2) ) ) 

CXc{2)=-(XP(?i**?/yc(i)-i,/>f '1>  **2-UP*LC AF  (2)/(SCFT(L0AP(2)**2 
1  +  LDAP  (  2 ) **2 > ) ) 

CXC(7)=-{-XF(2)*XP(3)/XF(l )-LF*LOAo (  7)  / <S CRT (LDAP ( 2) +  +  2 +  LD A F  (  3 ) 
i  * *2)  )  ) 

17  2  OF(l)  =-XE  { 2) 

CXF(4>  =-XE(3)/X£(i) 

CXF(l)  =-X»(2) 

C  X  P  ( 4 )  =-XP(3)/XP(l) 

CLCA"(1)=-((LCAE(2)  *XF,(3)**2)/yE(l)**2“2.*LCA£(2)/XE(l)**3 
1  -(L0A2(3)*XE(?)*XE  (3)  )/XE  <D#*2+L'3AE(4)  *XC  (3)/XE  (  1)  **2) 

CL  CAE (C)s-(-LCAE (1) +  LCAE  (3) «XE  <3)/XE(i> ) 

ClCAE(o)=«(-{2«*LOAE{2)*Xf(2))/XF(i)+LOAE  (3)*XE(2)/XE(i)-LQAE(4 
1  / XE ( i  )  ) 

CLCAF(l)=-(  ( L  C  A  P  (  2 )  *  X  P  (  3  )  *  *  2  )  /  X  P  ( 1 )  **?-£.  *  L  C  A  P  (  2 )  /  XP(1)  *  *  3 
.1  -(L0Mp(3)*Xc(2)*XP(3))/XF(n<Jf2+LJAD(4)«XF(3)/XC(l)»J»2) 
CLCAP(2)=-(-LCAP(l)  +LOAF(3)*XF(3)/XP(l) ) 

CLCAC (3)»- (-(2«*LD4  F  <2)*>F  (2  )  )/XP(i) +LDAF  (2) *XP (2) /XP(1) -LCAF  (A 

1  /X° ( 1 )  ) 

IF  (Mi E C  •  1 )  PC  TO  4  C 
CC  3n  1=1,4 
SXC(I) =XE(I) 

S X F ( I ) =  X°  (X) 

5CXE(I) =0XE(T) 

30  SCXP(I) =CXP(T) 

CC  31  11*1,3 
SLE(Il) sLCAE(Il) 

SLF(Il) =  L  CAP  ( II ) 

SCLOE(Il)  =DL0AE  ( ID 
7i  SCLr7D(Il)  =  CL0AP(I1) 

CC  32  12=1,4 

XE (I?) =  XE (12) +  OX  E ( I  2  )  *CT 

32  Xf (12) =XF (12) +GXF (I  2 )*C7 
CC  33  I ?=  1 , 3 

LCAE(I3)  =  LOAr(I3)+'OLCAE(TDJ»CT 

33  LCAP(I3)=LDAF  (13)  +DLCAP  ( 1 3  )  *  C  T 
►  *1 

CC  TO  1C1 
4n  CC  41  J  =  1  ,  4 

XEAVG(J)=(SrXE(J)+7XF(J>>/<. 

41  XFAVG(J)  =  (SDXF (  J  )  +  Q  X  F  (J)  )/?» 

CC  42  Jl=l,3 

LEA\/G(Jl)  =  (SCLOE(Ji  ) *01 C A  E  ( w li ) /?. 

42  LFAVG ( Jl)  =  (SuLOP (J1 ) 4TLCAF  (w  1  )) /2. 

CC  ?0  K  =1  ,4 

XE(K)=PXE(Ki+XEAVC-(K)*OT 

50  XF  (K) "SXP  (K)  +  XPA  VC ( K ) *CT 
CC  51  K  1  =  1 »  3 

LCAE(K1)=SLE ( K 1 ) +LEAVG(Kl) *  C  1 

51  LCAP(K1)=SLP(K1)+LPAVG(K1)*C1 
IF  (N.  EC.  ICO  GO  TO  fJ 

N  =  N  +  1 
R  =  ” 

GC  TO  lCl 
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63  Y  U)  =  m<l>-*P<1)> 

Y (?) = (X-  (2)  -XP( 2) ) 

Y  (3)  =  (  Yu  U)  -XP  (  3)  5 

Y  ((*)  =  (  Yr  U>  -XO(U)) 

FTS=SQFT(Y(1)**2  +  YC«<)**2) 

VFAT=Y  (2  )/Y  (2) 

V'fOT  =  3CPT  (Y  ( 2 ) *  *2  +  Y  C3>**2  > 

T  au=TAt+  0. 031 

VPrTC  (6,26)  Y (1) , Y  (2)  , Y ( 2)  ,Y (u) ,2TS, VR M  ,  VTCT,TAl 
26  FC3MAT  (1F8E15.6) 

IF  (TAU.CT. 3.02)  GC  TC  7C 

Y  -3 
N  =  t 

GC  TO  10 C 

73  IF  <FL2. IT. ?Eb. )  GC  TO  SC 
F  c£=FE  E- 5 , 

T  A  i  J  =  0  . 

XF ( 1 ) = i « 

XF(2) =C. 

XP(3)  =1, 

XP(4)  =C. 

GC  TO  13 
33  CONTINUE 

kFITE  (6,501) 

501  FC=?MAT  (*S*1 
E  NO 
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Linear  Integration  Program 

FSOf.^A1'  PAIN  (INPUT,  ClTd*  ,  TAP£6=I  N^L'T  ,  1  AF£e=0'JTPLT) 
tl^ENSIOK  Y  (1*)  ,rY<  4)  ,  SY  a  )  ,SCY<4)  ,  LOA  (  U  )  ,CL0A  (3>  ,  YAVG  <  4)  , 

1  L.JAvG(3),SLCA{?>,  SDLCAd) 

REAL  LCA,LPAVO 

SlAD  <5,  m)F^,Y2C,R,E7,T/!l 
1JS  FORMAT  <5rl2.G> 

UPIT^  (6  ,5C  .  ) 

501  F  C  "  A  T  (*T*> 

L  =  I 

63  IF  (L.E0.73)  GO  TO  66 

KV  =  1 

F  F 1  =  (F£c*3«  lsl“a)  /  16  0. 

TAUPHI  =STN  t  eM )  /CO  S  ( CHI ) 

LCA(l) =  SIN(FHI) 

LCA ( 2 ) =  C. 

LCA(3)=  C. 

LCA( 4)  =  COS  (PHI) 

RRITE  (6  ,  i C  ) 

10  FORMAT  (*  FHI  -  Y20  R  6  T 

1  7AU*> 

RFTT£  (6,ll)F£c,Y2C,R,FT,TAL 

11  FORMAT  (iP^-16.5) 

C^J.  0  J *3  01 

F=0 

N=1 

Y  (1)  =R*«IK(PH) 

Y  (2)  =  Y  2  0 

Y  (3)  =  P',SIN'(PHI)  -  Y2  1*TANFHI 


Y  (4)=R»CCS  (FHI) 

KFITE  (6,11;) 

F  CRH  A  T  {*  Y ( 1 ) 

Y  (2) 

Y  (?) 

YU) 

RTS 

Y  (  2  )  /  Y  ( 3 ) 

VTCT 

T  A  0  *  ) 

PTS=S0ST(Y(i;**2+Y(A)**2> 

\< p A T  =  Y  (2)/Y(3) 

VTCT  =  SCRT (Y (2) *  *  2  *  Y  ( 3  )  *  *  2  ) 

K-IT£  (6,20  Y  (i?*Y(2)  ,Y  (3)  ,Y(4),RTS,VFJT,V7CT,TAU 
2o  F  CPU  AT  ( 1 P6  i.  1.r  .  5 ) 

101  IF  (TAL.GT.  ’.)  CO  TC  1C2 

C Y ( 2 ) -  -2.*Y(o)-Y(l)+ET*STN  (FHI) 

CY(3) =  Y  (2) +  FT  *0  OS  (FHI) 

C-C  T0  lio 

112  CY(2)  =  -2»*VC*)-Y(  1)-PT*(LCA(2)/(SQRT  (L"A(2)**2+LCA(3)**2))) 

CY(3)  =  Y  ( 2 )  -  E  T  *  ( L  C  A  ( 3  )  /  <  S  C  F  T  (LCA(2)**2-*ICA  ( 3  )  *  *  2  )  )) 

103  CY( x)  =  -Y (25 

r*  v  /  i .  %  •»  .  v/  /  i  \  •  u  «  «  t 

L  •  \  t  -  ~  l  U‘M  l  Ul 

CLIA ( 1 )  =  LOA (2) -LH  A  (4) 

CL1A (25=  LPA (1) ~L0 A  {  3) 

CLQA ( 3  )  =  2. *LOn  (2)  +  LDA  (  U 
IF  (M.IC.l)  CO  TO  4[ 

CC  30  1=1,4 
SY(I) =  Y(I) 

30  SCY(I) =0Y(I) 
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%  *  /n  *  o 

IWW  V4.fi, '  fJA-i 


CC  3 i  11=1,3 
5LQ5(  ID  =  L  0  A  (  1 1) 

31  -  SCLOA  <  111=  3LDA  (ID 

CC  72  12=1,4 

32  Y  (1 2  >  =  Y  (12) +  OY  (12 ) *OT 
CC  ■*?  13=1,3 

33  L  D  A  ( 1 3  )  =  LDA(Il) +OLCA (13)  *C  I 
K  =  1 

C-C  TO  10? 

43  CC  41  2=1,4 

■+1  YA\/G  ( J  )  =  (SOY  (J)  +0Y  <  J)  )/?* 

CC  42  21=1,2 

42  LCDG<J1>=  (SOLOM  21) +CLCA  Ui)  )/2. 

CC  53  u2  =  l »  4 

53  Y(J2>=  5Y(  J2)  4-YAVG  (22)  ¥C7 

CC  51  23=1,2 

51  L  C  A  ( J  3  )  =  SLCA(J?)+trAVGU2)*CT 
IF  (N.En.iO-  GO  TO  eo 
MN+1 
F  =  Q 

GC  TO  132 

53  RT3=  SC^TIY  (i)**2+Y  (4)**2) 

V FAT  =  Y  (2)/Y (3) 

VTQT“  <ORT(Y  (2) **2  +  Y  (7)**?)- 
T  A'J  =  TAL  +  C.  C  jCI 

WRITE  (6  ,9  0)  Y(l)  ,Y<2)  ,Y<2)  ,Y<4),^TS,VHT,VTCT,TAU 
93  FORMAT  (lPer15.5) 

IF  MM.SC.2"')  GC  TO  70 
K  =  0 
N  =  1 

rM=MM+l 
GC  TO  101 
70  FC2  =  FC  6-5 • 

TA'J=  J . 

L  =  L  + 1 
GC  TO  63 
66  CCNTINIE 

WRITE  (6,501) 

5)1  FORMAT  DSD 
tKO 
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Linear  Closed-Form  Equation  Program 

FROGRAh  f-AIr,  (  INPUT  .OUTFIT  ,  T  AF  £5  =  1  NPl'T  ,  T AFE5 =CUT PIT) 

REAL  K1 , *2» <3, KU , h >  ,N 
READ  (E,  10  0)  FE'i’i  2 1. ,  R  ,  9  T  »RC»TAU 
130  F  CRM A  T  (EEin.5) 

L  =  1 

63  IF  U.E0.2)  GO  TO  6f 

F 1-1  =  lFLc*y.  1 59)  /13G. 

TANFHI  =  51K'  (FrI)/CC3(FFT) 

A  A  =  SIN  (PH T) 

EG  =  2.*SIU(FH7>*CCS<FH7) 

CC  =  -5.  /6.  +2.*COS  (FhT  )**  c  -»SIN  (PHI  )**2*CCS  (FHI)  **Z 
CC  =  CCS  (PHI) 

A  1  =  - A  A 

91  =  +CO 

Cl  =  - (CC/AA-5^*»2/  (  A. * AA** 3  )  >/2. 

Cl  =  -  (  <?.*FH*CS>  /  (  2  .*AA**3  >-  (7.*:H**2)  /  (fc.'AA**1;)  )/6. 

A  2  =  -  D  C 

92  =  -9c/  (2. *CO) 

C 2  =-(-CC/DC-P «»**?/  (4,*CC^^-3))/2, 

C  2  =  -((3.*r:r-’*3)/(E.*Cr««e)4t?.*0n*CC)/»2.J't:[:4*2n/E. 

K1  =-(2.4A2-Gl-4.4C2+£.*Cl)  ,ET 
K 2  =  -( 2. *??-?. *C1- 12. *D2) *ET 
K3  s-  (  2. *C2-2. *C1)  *ET 
K4  s-(2.*D2)*OT 
C  =  Y  2  T—  K  1 

F  =  - (  !, *R+~T) 4 AA+ 2 . *TANck  I*Yr0-K2 
MP  =  ( A. ^P+nT) *AA-2.*Y2.*TPhFHI+K2 

rv  =  -  (2.*P  +  ->n  »AAf  r  ANFFT#>2  c-K2 
S  =  R4CQ-?.*Y2  1+2. 4  K 1 
WRITE  (6,109) 

139  FORMAT  ( 4  nHI  Y23  9 

1  *> 

WRITE  (6,11'’)  FEE,Y2G,F,ET 
11C  FORMAT  (1FF-19.5) 

W^ITE  (6,129) 

129  F  CRHA  T  (*  Y1  Y2  Y3 

1  RV?  Y2/Y2  TACM 

K  =  1 

64  IF  (M.E0.21)  GO  TO  F5 
T  =  T  AU 

T 1  =  T  AU 

Y  1  =  P4  CO  S  (T  I)-04SIN  (Tl)-KlM-K?4T44?/2.-K?4T443/3.-K44T444/A 
1  +  MM 

Y  t  -  P  *  S I N  ill)  *  C  C  0  S  <  T  i  >  -,*!  +  »<.< 4  T  +  K34T4-  2+kW-7443 
Y2=-P*ruS(Tl)+n4STC.  (7i)  *  t  K  1  -t[T4A2)  4T  +  (‘''*5T4E2)4T442/2, 

1  +(K3  +  rT  4  C  2 )  4  T  4  *  3  /  -3  .  +  (kHFT*C?)*T**M./u.*h 
Y4-2,4C4SINni)42.*C*rrsni)-(2  4<1 4-0  T  4  A  2  )  4T442/2  .-(2  ,  *K2 
1  tOT* 2  2)  *T**7/-i.  -(  2  .  4K?4ET  4C2)  4T44  4/1  2  .  -  (2,*K'-  +  9T  *02) 

-  *  t  *  *  y  /  ^  : .  +  ( f-  m  -  ( ; )  *  t  +  s 

rT3=SOcT  (Yi4424Y44  4  2) 

RATY23-Y2/Y 3 
UTOr-^r^y  ( y  ?  *  *  2  4-  Y  3  ♦  4  2 ) 


jgWPPW»«WW‘*«1gr  nnrw  — —  -r^r km urnRov.ir:  -  »>- 

f?  '"  ■  -" 


f 
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^  K  X  T  E  {F.,13'')  Yi,Y2,Y3,Y‘(,R15-,r>ATY23,'r/ 
13  u  F  CRH  AT  (iP£u  i^.E) 

F-Mfi 

'AU-TAlf  C. u M 
CC  TO  ti* 

6«»  CCS'TTNIE 
T AU  =  C. 

FEE=Fc5-5. 

L=U  +  l 
CC  TO  e3 
&e>  CCi'J  T  INI? 

£  NO 


L  ,VTCT 
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